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Abstract 

An alternative derivation is presented of the infinitely many exceptional Wilson 
and Askey- Wilson polynomials, which were introduced by the present authors in 2009. 
Darboux-Crum transformations intertwining the discrete quantum mechanical systems 
of the original and the exceptional polynomials play an important role. Infinitely many 
continuous Hahn polynomials are derived in the same manner. The present method 
provides a simple proof of the shape invariance of these systems as in the corresponding 
cases of the exceptional Laguerre and Jacobi polynomials. 

1 Introduction 

In a previous paper [1] , we have derived infinitely many exceptional orthogonal polynomials 
related to the Wilson and Askey- Wilson polynomials [21 |3l H] as the solutions of infinitely 
many shape invariant [5j thus exactly solvable discrete quantum mechanics in one dimension 
[SI [TJ [HI [9l [To] . The corresponding Hamiltonians are the deformations of those for the Wilson 
and Askey- Wilson polynomials [6l [H [9] in terms of a degree i {£ = 1,2, . . .) eigenpolynomial 
with twisted parameters [1]. In discrete quantum mechanics, the Schrodinger equation is a 
second order difference equation, which reduces, in known cases of exactly solvable exam- 
ples, to the difference equation satisfied by the polynomials belonging to the Askey scheme 
of hypergeometric orthogonal polynomials and their g-analogues. As stressed in our previous 



publications, various concepts and formulas of these Askey scheme polynomials can be un- 
derstood and formulated in a unified fashion through the framework of quantum mechanics; 
that is the eigenvalue problem of a hermitian (self-adjoint) linear operator (the Hamiltonian) 
in a certain Hilbert space. The Hamiltonian corresponding to a particular polynomial, say 
the Wilson or the Askey- Wilson polynomial, is specified by a choice of the potential function 
(I2.9p as in the ordinary quantum mechanics. 

In this paper we present an alternative derivation of the infinitely many exceptional 
Wilson and Askey- Wilson poljTiomials in terms of Darboux-Crum transformations [Til 112] 
intertwining the Hamiltonians of the original Wilson/ Askey- Wilson polynomials with that of 
the corresponding exceptional polynomials. The same method is applied to the continuous 
Hahn polynomials [1] to construct the new infinitely many exceptional continuous Hahn 
polynomials indexed by a positive even integer £. This is a discrete quantum mechanics 
version of the recent work ^3] which derives the four families of infinitely many exceptional 
Laguerre/Jacobi polynomials [El [15] in terms of Darboux-Crum translations intertwining 
the Hamiltonians of the exceptional polynomials with those of the well known exactly solvable 
Hamiltonians of the radial oscillator and the Darboux-Poschl- Teller potentials [161 CZl [18] . 

The concept of exceptional {Xi) polynomials was introduced by Gomez-Ullate et al [19] in 
2008 as a new type of orthogonal polynomials satisfying second order differential equations. 
The Xi polynomials start with degree £ (£ = 1,2, . . .) instead of the degree zero constant 
term, thus avoiding the constraints by Bochner's theorem [20]. They constructed the Xi 
Laguerre and Jacobi polynomials [19], which are the first members of the infinitely many 
exceptional Laguerre and Jacobi polynomials introduced by the present authors [TH [TSj in 
2009. Quesne gave a quantum mechanical formulation based on shape invariance [2I| and 
later she introduced another set of X2 polynomials [22] . 

This paper is organised as follows. In section two we review the discrete quantum mechan- 
ical systems of the continuous Hahn, Wilson and Askey- Wilson polynomials in § 12. II together 
with those of the corresponding exceptional polynomials in § 12.21 The properties of the de- 
forming polynomial ^£ are discussed in some detail in § 12. 3[ They are important for deriving 
various results in the subsequent section. The main part of the paper, the Darboux-Crum 
transformations intertwining the original and the deformed systems are discussed in section 
three. The final section is for a summary and comments including the annihilation/creation 
operators and Rodrigues type formulas for the exceptional polynomials. 
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2 The original and deformed systems 

Here we first recapitulate the shape invariant, therefore solvable, systems whose eigenfunc- 
tions are described by the orthogonal polynomials; the continuous Hahn, Wilson and Askey- 
Wilson polynomials |4], to be abbreviated as cH, W and AW, respectively. See [6] and [9] for 
the discrete quantum mechanical treatment of these polynomials. It is based on the factori- 
sation method and generalisation of the Darboux-Crum transformations [HI [12]. See also 
the related papers [22] ■ The factorisation of the Askey- Wilson polynomial was also discussed 
in [21] . Then the deformed systems corresponding to the exceptional Askey type polynomial 
are summarised in § 12.21 

2.1 The original systems 

Here we summarise various properties of the original Hamiltonian systems to be compared 
with the deformed systems which will be presented in § 12. 2[ Let us start with the Hamilto- 
nians, Schrodinger equations and eigenfunctions (xi < x < p = —i-^): 




(2.1) 
(2.2) 
(2.3) 
(2.4) 



H{X) = ^(A)U(A), 

■H{X)Mx; A) = SnWMx; A) (n = 0, 1, 2, . . . 
(pnix; A) = (j)o{x; A)P„(?7(a;); A). 



The set of parameters A = (Ai, A2, . . .) are 



cH 



A = (01,02), ReOj > {i 



1,2) 



(2.5) 



W 



A = (ai , 02 , as , 04) , Re > (z = 1 , . . . , 4) , 
{a*, 03, 04} = {tti, 02, cis, 04} (as a set) 



(2.6) 



AW 



g-^ = (01,02,03,04), |ai| < 1, (i = l,...,4), < q < 1, 
{o*, Og, 03, O4} = {oi, 02, 03, 04} (as a set). 



(2.7) 



where ) = (^q-^^^q^'^^ . . .). The the sinusoidal coordinate r]{x) is. 




Xi = —00, ^2 = 00, 7 = 1 : cH 
Xi = 0, X2 = 00, 7 = 1 : W 
Xi = 0, X2 = 71, 7 = logg : AW 



(2.8) 
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The potential function V{x; A) and energy eigenvalue Sn{X) are 



v{x; a; 



def 



(ai + ix){a2 + ix) 

(2ix(2ix + l))~'nj=i(% + 
((l_e2-)(l-?e2"))-^nL(l 



tx) 



cH 
W 
AW 



def 



n{n + bi - 1), 
n{n + &!-!), 
(g-"-l)(l-64?"-i) 



def 



61 = ai + 02 + a* + al 
Oi = ai + a2 + 03 + 04 

, def 

04 = 01020304 



cH 
W 
AW 



(2.9) 



(2.10) 



Throughout this paper we consider the potential functions, eigenfunctions, etc as analytic 
functions of x in the complex region containing Xi < Rex < X2- We use the ^-operation 
on an analytic function * : / h-). /* in the following sense. If f{x) = Yl anX^, an G C, then 

n 

f*{x) =^ ^a*x", in which o* is the complex conjugation of o„. Obviously f**{x) = f{x) 

n 

and f{x)* = f*{x*). If a function satisfies /* = /, we call it a 'real' function, for it takes 
real values on the real line. 



The eigenfunctions are chosen real. 



and 



60 and P* = Pn- The main part 



def 



consists of an orthogonal polynomial Pn{v] ^ polynomial of degree n in rj: 

PniV; A) 

They are expressed in terms of the (basic) hypergeometric functions [1]: 



Pn{r]] 01,02,0*, 03) 
VF„(?7; 01,02,03,04) 
Pniv, 01,^2, 03,04!^) 



cH 
W 
AW 



(2.11) 



cH: p„(?7(x); ai,a2,a*,o;) 

def .„ (o-i + aj)„(ai + 02)0 



nl 



3-^2 



-n, n + Oi + 02 + o]; + O2 — 1, Oi + 2X 

Oi + Oi, Oi + O2 



1 , (2.12) 



W: iy„(r7(x); 01,02,03,04 

def 



(Oi + 02)„(ai + 03)„(ai + 04)n 

n, n + X]t=i o.j — 1, fli + ix, Oi — ix 



X 4F3 



Oi + 02, Oi + O3, Oi + O4 



AW: pn{r]{x)]ai,a2,a3,ai\q) 



def 



o-i ''(ai02, 0103, O1O4 ; q)n x 4^3 



q 01020304^" ^, OiC*"^, OiC ''^ 



(2.13) 



q;q], (2.14) 



Oifl2) Q-l^lS) O1O4 

in which (o)„ and (o; q)n are the Pochhammer symbol and its g-analogue. They are sym- 
metric in (oi, 02) for cH and in (oi, 02,03, 04) for W and AW. 
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These Hamiltonian systems are exactly solvable in the Schrodinger picture. Shape in- 
variance [5] is a sufficient condition for the exact solvability. Its relations involve one more 
positive constant k: 

^(A)^(A)^ = nA{X + (5)U(A + 6)+ ^i(A), (2.15) 

defHil) :cH dcffl :cH,W 

I (iiil) :W,AW' \ q-' :AW ' ^''^^^ 

or equivalently, 

V{x - q-, X)V*{x - ; A) = V{x; A + 6)V*{x -ir,X + S), (2.17) 

V{x + q; A) + V*{x- z|; A) = k{V{x; X + 6) + V*{x; X + 6)) - Si{X). (2.18) 

It is straightforward to verify these relations for the above potential functions for the cH, W 
and AW cases (12.91) . The groundstate wavefunction (f)o{x; A) is determined as a zero mode 
of the operator ^(A), A{X)(f)o{x; A) = 0, namely, 



'V*{x - q; A) - %l\ A) = ^V{x + z^; A) <l>^{x + A), (2.19) 
and its explicit forms are: 



T(ai + ix)V{a,2 + %x')V{a\ — %x^V{a\ — ix) : cH 



,^[x-X)={ J{V(2%x)V{-2ix)Y^X\]^^V{aj^ix)V{aj-%x) :W . (2.20) 



(e2" ; g)oo(e-2- ; q)^ n;=i(%e- ; qY^{aje-'- ; q)-^ : AW 
We introduce an auxiliary function V5(x) 



V9(x) =^ I 2x 



1 

2: 
2 sinx 



cH 

W , (2.21) 
AW 



which possesses the properties: 



(x; X + 5) = ip{x)^V{x + z|; A) 0o(x + z^; A), (2.22) 

-1 - ^ 7) - 



y(x;A + (5) = ^-^ ^^ _^^^^V (x-z|;A). (2.23) 



The action of the operators ^(A) and .4(A) ^ on the eigenfunctions is 

A{XY^{x- A) = /„(A)0„_i (x; A + (5) , (2.24) 



A{Xy(l)n-i{x;\+ S) = 6„_i(A)0„(a;; A). 



(2.25) 



The factors of the energy eigenvalue, /n(A) and 6„_i(A), £^„(A) = /n(A)6„_i(A), are given 

by 



/n(A) 



dcf 



n + bi-1 
—n{n + 61 — 1) 

5t(5-"-l)(l-645"-l) 



cH 

W , 6„_i(A) 
AW 



dcf 



n 

-1 

n 

q 2 



cH 
W . 
AW 



(2.26) 



The forward and backward shift operators J^(A) and B{\) are defined by 



J^(A) = 0o(^; A + 5)-' o ^(A) o <^o(x; A) = i^{x)-\eiP - e^i^) 



(2.27) 



^(A) = (/)o(x; A)-^ o A{\)^ o 0o(a;; A + 6) = -i{V{x] A)e^f - y*(x; A)e-^f)(^(x), (2.28) 



and their action on the polynomials is 



T{\)Pn{r]{x)- A) = U\)Pn-i{r]{x)\ A + 5), 
B{X)Pn-i{r]{xy, X + 8) = bn-i{\)Pn{v{x)] A). 



(2.29) 
(2.30) 



The second order difference operator ?^(A) acting on the polynomial eigenfunctions is defined 

by 



n{X) = B(A) J-(A) = (/.o(x; A)-^ o 7^(A) o (f)o{x] A) 

= V{x; X){e^P - 1) + V*{x; X){e-^P - 1), 
n{X)P„{rjixy, A) = £„(A)P„(r7(x); A). 

In conventional terms, this is the difference equation determining the polynomials: 



(2.31) 
(2.32) 



V{x; A) {Pnivix - 27); A) - Pnivix); A)) + V*{x; A) {Pn{vix + ^7); A) - Pnivix); A)) 
= £„(A)P,(?7(x);A). (2.33) 



The orthogonality relation is 
/ (j)o{x; A)^ -Pn(^(2;); X)Pm(r){x); X)dx = /i„(A)(5„ 



hn{X) 



dcf 



27rn?j=i r(n + ai + a*) ■ (n!(2n + fei - l)r(n + b, - 1))"' 
27rn! (n + 61 - 1)„ ni<i<i<4 r(?^ + 0. + %) ■ r(2n + 61)-^ 

27r(64g"-^ Q)n(&4Q'^ g")oo(?+^ <Z)oo' ni<i<,<4(««%5"; ?)^' 



cH 
W 
AW 



(2.34) 
(2.35) 
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As shown in detail in [8] these Hamihonian systems are exactly solvable in the Heisenberg 
picture, too. The positive/negative frequency parts of the exact Heisenberg operator solution 
of the sinusoidal coordinate r](x) provide the annihilation/creation operators a*^^-'(A) which 
map the eigenfunctions to the neighbouring levels 

a'^^\X)(f)n{x; A) oc (l)n±i{x; A). (2.36) 

This is a disguise of the three term recurrence relations of these orthogonal polynomials. The 
above relations are to be contrasted with the actions of the operators A{X) and A{X — Sy 
( I2.24p -( l2.25p . which maps <pn{x; A) to the neighbouring levels of shifted parameters \± d. 

2.2 The deformed systems 

Here we recapitulate the Hamiltonian systems of the exceptional Wilson and Askey- Wilson 
polynomials, which were derived by the present authors [1] in 2009. The exceptional contin- 
uous Hahn polynomials are new. For each £ = 1, 2, . . ., we can construct a shape invariant 
system by deforming the original system {i = 0) in terms of a degree i eigenpolynomial ^eif]) 
of twisted parameters. We restrict the original parameter ranges (I2.5l) - fl2.7p as follows: 

cH : Qi > 0, £ : even, (2.37) 
W: ai,a2GM, {0.3, 04} = {0,3, ^4} (as a set), 

0<aj<Reak (j = 1, 2; A; = 3, 4), (2.38) 
AW: ai,a2GM, {^3, O4} = {as, (as a set), 

l>aj>\ak\ (j = 1,2; A; = 3,4). (2.39) 

In terms of the twist operation t acting on the set of parameters, 

'(^)^^'{(:a;;-1a3.a.) iw.Aw 

the deforming polynomial C.iirj) is defined from the eigenpolynomial Peir]): 

Uv;^) = Pi{vA{^+ii-l)S)). (2.41) 

We may need to restrict parameters further in order that ^e{r]{x); \) has no zero in the 
rectangular domain xi < Rex < X2, \lmx\ < which is necessary for the hermiticity of 
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the Hamiltonian. The potential function, the Hamiltonian and the Schrodinger equation of 
the deformed system are: 

A) = A + M) ^^(,(, _ ,.)^ ^) 6(.(x);A + 5) ' ^^'^^^ 

V;i.; A) ^ V^i.; A + M) I a! ^^^J^ M ^^^^^ ' (^-^^^ 

+ A) ^^(/^(a;); A + d) 



^,(A) = t{e"rPJv;{x; A) - e-^^'yi^(^). 



^,(A)t -z(v/l^(^e?^ - ^V^/(x; A) e-i^'), (2.44) 
H,(A) =^,(A)U,(A), (2.45) 
niiX)<Pi,nix;X)=£i,„iX)Mx;X) (n = 0, 1, 2, . . .), £,,„(A) = £„(A + M). (2.46) 

The continuous Hahn polynomials are real polynomials defined on the entire real line. There- 
fore the odd degree members have at least one real zero, whichever coefficients we may choose. 
This is the reason why £ is restricted to even integers in the cH case fl2.37p . 
This system is shape invariant: 

AeWMXy = kMX + d^MX + S)+ f,,i(A), (2.47) 

or equivalently, 

Ve{x - X)V;{x - A) = Vfix; X + 5)V;{x - 27; A + 6), (2.48) 

V,{x + if, A) + V;{x - if, A) = K{Ve{x; X + d) + V;{x; X + d))- ^,,i(A). (2.49) 

Proof is straightforward by direct calculation. In order to derive eq. f l2.49p . use is made of 
the two properties of the deforming polynomial fl2.68l) - fl2.69l) and the factorisation of the 
potential f l2.70p . For another simple proof of shape invariance, see the discussion in the final 
section. The eigenfunctions are 

Mx; A) Mx-A + iS) ^2.50) 

Vdivix + if; X)Ce{v{x - if; A) 

<j),,nix; A) = X)P,,n{vix); A) (n = 0, 1, 2, . . .). (2.51) 

Here Pf,„(?7; A) is a degree i + n polynomial in 77 but P£ n{f]{x); A) has only n zeros in the 
domain xi < x < X2. The explicit forms of Pe^niv) were given by eqs. (42)-(44) in [Ij, with 
eqs. (66)- (68) for W and eqs. (80)-(82) for AW. Here we present much simpler looking forms 
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of them, which encompasses the new cH case, too: 

A) = . ~' ^ ^ U{x; X + i5)Uv{x + ^i); A)F„(r7(x - ); A + M + 
/^,„(A)^(x) ^ 

- vlix- A + iS)Uvix - ^i); A)P„(?7(x + qy, X + iS + 6)), (2.52) 



where 1^1(2;; A), fi^nW and 6 will be defined in (12. 71 p . (I2.72p and (I3.12p . This is one of the 
main results of the present paper which is derived in §[3l Its lowest degree member is the 
degree i deforming polynomial itself of the shifted parameters 

Pe,o{v{xy,X)=Uv{xyA + S), (2.53) 

which is obtained by (I2.68p . It is straightforward to verify that the groundstate eigenfunction 

0, o(x; A) = M-;^ + iS)Uv{^yX + 6) ^2.54) 
V^eivix + i^)]X)^£{'n{x~i^)]X) 

is the zero mode of the operator ^^(A), Ae{X)(f)ifi{x; A) = 0. 
The action of ^^(A) and ^^(A)''" on the eigenfunctions is 

Ae{X)(f)i^n{x; A) = /£,„(A)0£,„_i (x; X + d), (2.55) 
Ae{Xy(j)e,n-i{x] A + <5) = bi^n~iW<Pe,ni^; A), (2.56) 
knW = /n(A + M), b,,n-i{X) = 6._i(A + iS) . (2.57) 

Like the corresponding formulas of the original systems (I2.24p - (l2.25p . these are simple con- 
sequences of the shape invariance and the normalisation of the eigenfunctions. In the next 
section, we will derive these formulas through the intertwining relations and without assum- 
ing shape invariance. 

The forward shift operator J^e{X) and the backward shift operator Be{X) are defined in 
a similar way as before 

J-,(A) = M^-, A + S)-' o ^,(A) o ij,{x- A) 

= ^{xMvixyx) (^^(^^^ + ^ + - ^^(^^^ - ^ + ^)'"'') ' ^^-^^^ 

Bi{X) = Mx- A)-i o A,{Xy o ^,(x; X + 8) 

^ e.(.(xrA+^) K^^+^^^^-^^^^+^-^^^^^'^ 

- V*{x- X + m,{r^{x - iT)- A)e-^p)(/.(x), (2.59) 



and their action on the polynomial Pe,n{Vj ^) is 

^£(A)P,,„(r7(x); A) = /,,„(A)P,,„_i (?7(x); A + 5) , (2.60) 
-B,(A)P£,„-i(r7(x);A + (5) =6,,„_i(A)P,,„(?7(x);A). (2.61) 

The second order difference operator 'H^(A) acting on the polynomial eigenfunctions is de- 
fined by 

n,{X) = B,{\)F,{X) = ^,{x- A)-i o Ut{\) o ^,{x- A) 

^^4?^(a; + i|); A) ( _ ii[ri{x - %^)\\^ 6) 



i,{ri{x - e.(r/(x);A + 5) 

x6(^(a^ -^i); A) / ie,{'r]{x^ii)\\ + 6\\ 

+ V\x, A + ib r y) ^-' x - V r r /xVa^ ' ^2.62) 

^4r^(x + 2^); A) V 6(^(a^); A + d) / 

•H,(A)P,,n(r7(x); A) = ^,,„(A)P,,„(r7(x); A). (2.63) 



Again it is trivial to verify that the lowest degree polynomial Pi/i^{ri{x)\ A) = A + 5) 

is the zero mode of 'K(,{X): 

ne{XMv{xy,X + d)=0. (2.64) 

The orthogonality relation is 



^2 

2 



xi 



ifjiix; A) P^,„(?7(x); A)P^,m(?7(x); A)(ix = /i£,n(A)5„m, (2.65) 

' (2ai + n + £)(a2 + a^ + n + 2£-l) 



/i£,„(A) =^ /i„(A + M) X < 



(2ai + n){a2 + a*2 + n + i - 1) 
(ai + aa + n + £)(a3 + 04 + n + 2£ - 1) 



(ai + a2 + n){as + + n + £ - 1) 
_, (l-aia2g"+0(l-a3a4g"+^'~^) . 
^ (l-aia2g")(l-a3a4g"+^-i) 

2.3 Properties of the deforming polynomial 



cH 

W . (2.66) 



Here we present three formulas of the deforming polynomial ^^(77; A) fl2.67p -f B.69p . which 
will play important roles in the derivation of various results in section three, in particular, 
the fundamental results of this paper f l3.10p and (13. lip : 

(V{x- t(A + (£ - 1)6)) (e^^ -1) + V* {x; t(A + {i - 1)6)) {e'^^ - l))Uvix); A) 

= Se{t{X))Uv{^yA), (2.67) 
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^){vl{x- X + id)e'^P - v,{x- A + id)e-iP)Uvixy, A) = f.^oWUvix); A + S), (2.68) 



= bi,oW^e{'n{x)]X), 

where A), V2{x] A) are the factors of the potential function V{x] A): 

vi{x; X)v2{x; A) 



A) = -v^- 



a; 



i{ai + zx) 



cH 
W 
AW 



W2(x; a; 



def 



z(a2 + zx) 



-"n;=i(i 

The constants fe^nW and 6£,„(A) are given by 
2ai + n 



n 



i=3V 



cH 
W 
AW 



(2.69) 

(2.70) 
(2.71) 



dcf 



01+^2+'^ 

-qr- V(l - aia2q") 



cH 

W , 6,,„(A) 
AW 



dcf 



02 + + n + 2^ - 1 
as + + n + 2i — 1 



cH 

W .(2.72) 
AW 



The first equation (I2.67P is the difference equation for the deforming polynomial, which 
corresponds to (12.33^ . The eqs. ( I2.68p -( !2.69p are identities relating ^e{Vi ^) ^^iv'i A + (5). 
We remark that these two equations ( ]2.68p -( l2.69p imply ( 12.67p . In similar problems in 
ordinary quantum mechanics, the exceptional Laguerre and Jacobi polynomials, analogous 
identities play important roles in proving shape invariance and other relations pTf EHJ [13] . 
As shown in fl2.12p -( l2.14p . the continuous Hahn, Wilson and Askey- Wilson polynomials are 
expressed in terms of the (basic) hypergeometric functions 3F2, 4F3 and 403, respectively [1]. 
The identities f l2.68p - fl2.69p are reduced to the following identities satisfied by the (basic) 
hypergeometric functions: (a, a', ai, . . . ,04 : generic parameters) 



-a, a', ai + 1 + ix 

as, tti + ^2 + 1 
—a, a', «! + ix 
as, «! + a2 



ai + ix) 3F2 
(ai + 02) 3-^2 

«! + ix){a2 + ix) 4F3 
— (cti — ix){a2 — ix) 4F3 
01 + 02) 4F3 



1 ) + ("2 - ix) 3F2 



—a, a ,ai + ix 

as, Q!i + ^2 + 1 



(2.73) 



—I 
2^ 



—a, a', «! + 1 + ix, ai — ix 
ai + a2 + 1, «! + as, a\ + a4 
—a, a' , ai + ix,ai + l — ix 

«! + ^2 + li Ol + «3) + 



+ ix, ai 



IX 



ai + ^2, ai + as, ai + a^ 



(2.74) 



—I 



2sinx 



e-*"(l-aie'")(l-a2e") 4<; 



a ^ , a' , aiqe^^ , aie '' 
aia2q, diCts, aia^ 



q;q 
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- e''^(l - aie - age *'^) 4<^3 g; g 

q;q). (2.75) 

aia2, tti«3, ^ 

These identities can be easily verified by the series definition of the (basic) hypergeometric 
functions. 

3 Intertwining relations 

Here we demonstrate that the Hamiltonian system of the original polynomials reviewed in 
§ 12.11 and the deformation summarised in § 12.21 are intertwined by a discrete version of the 
Darboux-Crum transformation. This provides simple expressions of the eigenfunctions of the 
deformed systems fl2.52p in terms of those of the original system, which is exactly solvable. 
It also delivers a simple proof of the shape invariance of the deformed system. 

3.1 General setting 

For well-defined operators ^^(A) and ^^(A)''", let us define a pair of Hamiltonians 7i^^^(A) 

hJ+Ha) = i,(A)U,(A), h[-\X) =AeWM^)\ (3-1) 

and consider their the Schrodinger equations, that is, the eigenvalue problems: 

?if)(A)<^fj(x;A) =4^(A)0{5(x;A) (n = 0, 1, 2, . . .). (3.2) 

By definition, all the eigenfunctions must be square integrable. Obviously the pair of Hamil- 
tonians are intertwined: 

?i(+)(A)^(A)t = i,(A)ti,(A)^(A)t = MXYn[-\X), (3.3) 

i,(A)?if Ha) = AiWM^yM^) = n\'\x)M^)- (3.4) 

If Ae{X)(})^f^j^{x; X) 7^ and Ae{Xy (t)[J {x; X) 7^ 0, then the two systems are exactly iso- 
spectral and there is one-to-one correspondence between the eigenfunctions: 

4^(A)=4;Ha), (3.5) 
0(-J(x; A) (X AeW&{x; A), (x; A) oc i,(A)t0y (x; A). (3.6) 
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It should be stressed in the ordinary setting of Crum's theorem, the zero mode of ^^(A) is 
the groundstate of 'H'i^\X). In that case, V}t\X) and H['\X) are iso-spectral except for 
the groundstate (essentially iso-spectral) of ?i^^'*(A), since it is annihilated by ^^(A) [25] . 

In the following we will present the explicit forms of the operators ^^(A) and ^^(A)^, 
which intertwine the original systems in § 12.11 and the deformed systems in § 12.21 

3.2 Intertwining the original and the deformed systems of the 
polynomials 

The potential function Vi is the original potential function V with twisted parameters and 
multiplicatively deformed by the deforming polynomial ^f. 

U^- A) V{x- t(A + (£ - m f^f^;l^l,^\ (3.7) 
Vnx- A) = V* {x- t(A + (^ - l)^)) ^'^^/(^^(t)?A)^^ ' (3-8) 

dcf ./ 2r 



A,{\) = z(e^V\/;(x; A) - e-^^JVe{x- A) ) 



A(A)t "=^' -i{^V,{x- A) ei^ - ^JVl{x■ A) e~i^) . (3.9) 

It is illuminating to compare these potential functions f l3.7p -( l3l8ll with those of the original 
(12. 9 p and deformed (I2.42p -( l2.43p systems. Again it is obvious that the overall normalisation 
of the deforming polynomial ^£ is immaterial for the deformation. For this choice of ^^(A) 
and ^^(A)''", one of the pair of Hamiltonians 'H^^\\) (13. ip becomes proportional to the 
original Hamiltonian 'H(A) (12. 2p with X X + id + 5 and the partner Hamiltonian ""(A) 
is proportional to the deformed Hamiltonian l-Lt{X) (12.450 : 

nf\X) = ke{X){n{X + iS + S) + /,,o(A)6,,o(A)), (3.10) 
n'f\X) = k,{X){U,{X) + /,,o(A)6,,o(A)), (3.11) 

where ki{X) and 5 are given by 



. .def/l :cH,W odef r (2,-2) :cH 

^'^^^-\{a,a,q')-^ -.AW ' " \ (i, i, -i, -1) : W, AW ^^'^^^ 

The multiplicative and additive constants are common to the pair of Hamiltonians (I3.10p - 
(13. lip . These are the main results of this paper. Like the corresponding formulas in ordinary 
quantum mechanics [13], these fundamental results can be obtained by explicit calculation. 
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The three formulas in § 12.31 are essential. That is, the difference equation satisfied by the 
deforming polynomial ^f {rj{x) \ A), f l2.67p and the two identities relating the deforming poly- 
nomial ii{r]{x)] A) to its shifted one ii{r]{x)] A + <5), fl2.68p and (12.691) . 

It is instructive to verify that the zero modes of .A^(A) and ^^(A)^ do not belong to the 
Hilbert space of the eigenf unctions. In fact, the zero mode ^^(A) is 

A,{X)x = 0, x = Uvi^y. A)0o(2;; t(A + (£ - 1)8). (3.13) 

It has at least one pole in the rectangular domain xi < Rex < X2, |Imx| < ||7|, therefore 
it cannot belong to the Hilbert space. The zero mode of ^^(A)^ is 

- ii); A)&(,(x + ii); A) V(i- - if, t(A + l)S) 

which is obviously non-square integrable. This situation is the discrete analogue of the 
'broken susy' case in ordinary quantum mechanics in the terminology of supersymmetric 
quantum mechanics [T71 125] . 

Based on the results (I3.10p -( l3.1ip . we have 

4>^^^{x-\) = 4>n{x-X + l5 + 6), 0y(x;A) = (^,,„(a;;A), (3.15) 
^tnW = f^iW{£n{\ + M + ^) + /,,o(A)6,,o(A)) = k,{\){£,^^{\) + /,,o(A)6,,o(A)). (3.16) 

Then it is trivial to verify Ae,{\)(i))^^{x] A) 7^ and A(,{X)'^(j)'f^{x] A) 7^ 0. For, if one of the 
eigenfunction is annihilated by ^^(A) (^^(A)"!'), the left hand side of f l3.10p (( p.lip ) vanishes, 
whereas the right hand side is ki{\){£n{\ + C-5 + d) + o(A)6^^o('^)) times the eigenfunction, 
which is obviously non-vanishing. Note that £^„(A + 18 + 8) = £^„(A + 18). 
The correspondence of the pair of eigenf unctions (p^tni^) expressed as 

a/kHA) /^,„(A) ^JK|,[\)hl^nW 

We introduce operators Pi{X) and Bi{X) defined by 

P,{\) U^- X)-' o o </,o(x; X + i8 + 8), (3.18) 

4(A) 0o(x; X + e8 + 8)-' o 4= o M^; A). (3.19) 

'Ke{X) 
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The operators ^^(A) and Be{\) are expressed explicitly by using the concrete forms of 
V{x] A), ipe{x] A) and (l)o{x; A): 

fi{X) = ^ {v,{x; A + id)Uvi^ + ^i); A)ei^ - vl{x; A + id^vi^ - ^i); A)e-ip) , (3.20) 
^^(^) = X, 7\ A + - l)5)ei^ - v;{x; X+{e- l)S)e-iA . (3.21) 

Here is a technical remark on (13.201) . In deriving the explicit form of the operator ^^(A) 
in (I3.20p . one extracts the factorised potential function Vi{x) and vl{x) from the corre- 
sponding expression of the square root of the twisted potential ^/V{x; t(A + {i — 1)<5)) and 
^V*{x; t(A+ (^- 1)<5)). Thus the choice of the argument is a subtle problem, in particular, 
for the cH case, in which only one factor undergoes the sign change. In other (W and AW) 
cases, sign change occur in two factors and thus the effect cancels out. Let us consider the 
twisted potential of cH case 

V{x; t(A + i6)) = (-ai - ^+ix){a2 + ^ + ix). 

For positive x close to the origin (|x| <^ 1), we choose —ai — | + to have an argument 



close to +7T. Then its ^-operation —ai — | — = — (ai + | + ix) has an argument close to 



—TV. This would mean y — (ai + | + ix)"^ = —iio-i + | + ix) = ~Vi{x; A + id), instead of the 
naively obtained i{ai + | + zx) = fi(x; A + id). 

The operators -^^(A) and i3^(A) act as the forward and backward shift operators connect- 
ing the original polynomials Pn{v) ^-^id the exceptional polynomials „(//): 

f,iX)Pnivix); \ + i6 + 6) = fi,nWPM^y, A), (3.22) 
BeWPeAvi^y^ X) = b£,„(A)P„(r/(x); X + id + 6). (3.23) 

The former relation fl3.22p with the explicit form of ^^(A) fl3.20p provides the new explicit 
expression fl2.52p of the exceptional orthogonal polynomials, which is one of the main results 
of this paper. 

Other simple consequences of these relations are 

4? (A) = /c,(A)/,,„(A)b,,„(A), £rXX + id) = knWknW - koWk,o{X). (3.24) 

The normalisation constant /i^^„(A) (12.660 of the exceptional polynomials is related to 
that of the original polynomial /i„(A) (12.350 : 

heniX) = h^KiX + iS+~S) = ^5n(A)/o,n(A + M) 

' An(A) ^ faX)bo,n{X + iS) ^ ^ ^ 
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In the second equality we have used the exphcit forms of /i„(A) (I2.35p . Eq. fl3.25p is shown 
in the following way: 



Xl 



X2 



Xl 

X2 



(£) ^(+) 

(iv) 



dx Ae{\)(l)nix; X + i6 + 6)- Ai{X)(f).mix; X + i6 + 6) 

rx2 „ „ _ 

/ dx Ai{XyAe{X)(f)nix; X + iS + S)- X + i6 + 6) 

J X\ 

^It^W I ' dx X + iS + S)<Pm{x; X + £6 + S) 

J Xl 

k,{X)f,,n{X)h,n{X)K{X + 15 + 5)5nm. (3.26) 



Here we have used (IXTTj) and (IXm in (i), ([S2D and in (iii), flX^ and fICTD in (iv). 

In order to show (ii), we need to shift the integration contour to the imaginary direction. It 
is allowed if Vi{x\ X)(f)Q{x\ X + 15 + 6Y has no pole in the rectangular domain xi < Rex < X2, 
g ^ ^ |_ This condition is fulfilled if the deforming polynomial C,i{'r]{x); X) has no zero 
in the rectangular domain xi < Hex < X2, |Imx| < ||7|, which is indeed the case. 

3.3 Other intertwining relations 

It is interesting to note that the operator Ae{X) intertwines those of the original and deformed 
systems A{X) and Ae{X): 

Ai{X + S)A{X + iS + S)= Ae{X)Ai{X) , (3.27) 
^(A)^(A + iS + Sy = AiiXyAiiX + S). (3.28) 

In terms of the definitions of the forward shift operators J^{X) fl2.27p . J^e{X) (12.581) . -F^(A) 
fl3.18p . and B{X) (12.280 . Bi{X) (12.590 . the above relations are rewritten as: 

y/kfXX + S)M>< + ^)^(A + £6 + 6) = ^/^X}J^fXX)M^), (3.29) 
^/ ke{X) :Fe{X)B{X + i6 + 6) = ^/k^{X + 5) Bi{X)Pi{X + 5). (3.30) 

These relations can be proven by explicit calculation with the help of the three formulas of 
the deforming polynomial ^^(r/; A) (I2.67p -f l2.69p in § 12.31 

By applying ^^(A + 6) to ^TIM and ^^(A) to ^TM together with the use of ([32ID> 
fl3:28|) and fIXTTp . we obtain 



At{X)Mx-. X) = W ^^^^^ /n(A + £6 + ^) /An-i(A + ^) ^ x + 6) 

l^iW fi,nW 
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= fn{X + i6)^i,n-l{x-X + d), (3.31) 

A(A) V.,n-i(x; X + S) = J . ';[^^\. 6„_i(A + £5 + 5) . ^^'"^^^ c/>,.„(x; A) 

= &„_i(A + £(5)0,,„(x;A + (5). (3.32) 

In the calculation use is made of the explicit forms of k^(A), /^^„(A), /n(A) and &n(A) in 
the second equalities. This provides a proof of fl2.55l) - fl2.57p without recourse to the shape 
invariance. Likewise the above intertwining relations of the forward-backward shift operators 
fl3.29p - p.30p give the simple proof of f l2.60p - p.6ip . respectively, again without recourse to 
the shape invariance. 

4 Summary and Comments 

The Darboux-Crum transformations intertwining the Hamiltonians of the continuous Hahn, 
Wilson and Askey- Wilson polynomials with those of the corresponding exceptional polyno- 
mials are constructed in a unified fashion. This gives a much simpler expressions f l2.52p of 
the exceptional Wilson and Askey- Wilson polynomials than those given in a previous paper 
[1]. The exceptional continuous Hahn polynomials are new. See the recent work of Gomez- 
UUate et al [26]. It provides the iso-spectral Darboux-Crum transformations intertwining 
the Hamiltonians of the Laguerre polynomials to that of the radial oscillator and shows 
the shape invariance. 

The present paper is a discrete version of the recent work [13j which provides the Darboux- 
Crum transformations intertwining the Hamiltonians of the radial oscillator/Darboux-Poshl- 
Teller potential (the Laguerre and Jacobi polynomials) with those of the corresponding 
exceptional polynomials [HI [TSl EZ]- This offers a simple proof of the shape invariance of 
the i-th exceptional polynomials through the established shape invariance of the original 
polynomials as depicted in the following commutative diagram: 
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Ai{x] X + 6) 




established 
shape invariance 



shape 
invariance 



(xn{X + id + d) + c{X) 



Ai{x; X) 




original polynomial exceptional polynomial 

Two ways of proving shape invariance of the exceptional polynomials system. 

The Darboux-Crum transformation also supplies the annihilation/creation operators for 
the exceptional polynomial systems through those a*^^-'(A) for the original system fl2.36p 



(4.1) 



which map (j)e,n{x; X) to (l)e,n±i{x; A). The analogous formulas work for the exceptional La- 
guerre and Jacobi polynomials [13]. As shown in [HI [9] , the annihilation/creation operators 
together with the Hamiltonian constitute dynamical symmetry algebra of an exactly solvable 
system, for example, the g-oscillator algebra [29]. It would be an interesting challenge to 
clarify the structure of the dynamical symmetry algebras associated with the exceptional 
Askey type polynomials. Likewise the three term recurrence relations for the original poly- 
nomials are mapped to those of the exceptional polynomials [2H] . However, their significance 
and utility dbTG clS yet unclear. 

In [27], the shape invariance of the infinitely many potentials for the exceptional La- 
guerre and Jacobi polynomials are attributed to as many cubic identities among the original 
polynomials. The corresponding identities for the cH, W and AW polynomials are given by 
f l2.49p . These identities are quartic in and can be proven by using fl2.68p and f l2.69p . 

Rodrigues type formulas for the exceptional polynomials „ are obtained by multiple 
applications of the backward shift operators for the exceptional polynomials ^^(A) f l2.6ip 
or multiple applications of the backward shift operators for the original polynomials B{X) 
f l2.30p followed by the intertwining forward shift operator ^^(A) f l3.22p 



n-l 



P,,„(77(x);A) = n 



Be{X + kS) 



k=0 



'n-l-k 



'X + {l + k)5) 



Unix)] X+{n + 1)5) 



(4.2) 
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/,,„(A)f-i&n-i-.(A+(£ + A:)5) ' 

where 11^=0 '^k = clqCLi ■ ■ ■ fln-i is the ordered product notation of operators. In this con- 
nection, it is interesting to compare the forward-backward shift operators for the original 
poljTiomials f l2.27p -f l2.28p . the exceptional polynomials f l2.58p -f l2.59p and the intertwining 
ones fl3.20p -f l3.2ip . For the original polynomials, the forward shift operator J-'(A) fl2.27p is 
trivial, whereas the potential dependence is contained in the backward shift operator B{X) 
f l2.28p . For the exceptional polynomials, the forward shift operator ^-'^(A) fl2.58p is solely de- 
termined by the deforming polynomial and the potential function enters in the backward 
shift operator Be{X) fl2.59p . The intertwining ones are really twisted. The twisted part vi{x) 
of the potential function enters in ^^(A) fl3.20p . whereas the untwisted part V2{x) remains in 

Certain generating functions of the exceptional polynomials are easily constructed from 
those of the original polynomials through the main result f l2.52p . Suppose a generating 
function of the original orthogonal polynomials Pn is given by 

oo 

G{t, x;X) = J2 an{X)Pn{v{xy, A)^, (4.4) 

n=0 

in which an{X) is a constant. For explicit forms, see [1], eqs. (1.4.11)-(1.4.13) for cH, (1.1.12)- 
(1.1.15) for W and (3.1.13)-(3.1.15) for AW. Then eq. fl2.52p gives the corresponding gener- 
ating function of the exceptional orthogonal polynomials Pi^n, 

oo 

«n(A + i6 + ^) A„(A)P,,„(r/(x); A)^ 

(^vi{x; X + iS)Uv{x + qy, X)G{t, x-if,X + i6 + 6) 
- vlix; X + id)Uvix - qy A)G'(t, x + q;X + iS + S)^ . (4.5) 



n=0 

—i 

if{x 



It is well known that the Laguerre polynomials are obtained from the Jacobi polynomials 
in a certain limit. Likewise the Wilson polynomials are produced by a certain limit from 
the Askey- Wilson polynomials. The corresponding limiting relations for the exceptional La- 
guerre and Wilson polynomials are discussed in [15] and [28] , respectively. Here we comment 
on the limiting relations among the exceptional continuous Hahn polynomials and the excep- 
tional Wilson polynomials. The continuous Hahn polynomials are obtained from the Wilson 
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polynomials in the following limit: 

lim (— 2L)~"W„((x + L)^; ai — iL, a^ + iL, a2 — iL, a4 + iL) = n\ Pn{x; ai, 02, 03, 04). (4.6) 

L— s>oo 

By taking = + L and = {af - iL,af + iL,af - iL,af* + iL), and after 
appropriate overall rescaling, various quantities of the exceptional Wilson systems reduce to 
those of the exceptional continuous Hahn system in this L — )■ 00 limit. 

The idea of the infinitely many exceptional Laguerre and Jacobi polynomials [H] was ob- 
tained while studying various possibility of generating exactly solvable quantum mechanical 
systems from the known shape invariant ones, e.g. the radial oscillator and DPT potentials. 
Adler's modification |30l [31] of Crum's theorem is the most comprehensive way to generate 
infinite variety of exactly solvable systems from a known ones [321 125] . After the formulation 
of the discrete quantum mechanics version of Crum's theorem [331 1311135], its modification a 
la Adler is now published [36]. Its Appendix has many formulas reminiscent of those given 
in this paper. 

There are two types of discrete quantum mechanics. In one of them, as discussed in 
this paper, difference operators cause shifts in the pure imaginary direction [9]. The for- 
mulation of the other type of discrete quantum mechanics, in which difference operators 
cause real shifts, is provided in [10] . The corresponding eigenfunctions of the known solvable 
systems consist of the orthogonal polynomials of a discrete variable, for example, the (q-) 
Racah polynomials, [H [37] . It is a good challenge to construct the exceptional polynomials 
corresponding to these orthogonal polynomials of a discrete variable. 
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